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ABSTRACT 

The inner Galactic Bulge has, until recently, been avoided in chemical evolution studies due to extreme extinction and stellar crowding. 
Large, near-IR spectroscopic surveys, such as APOGEE, allow for the first time the measurement of metallicities in the inner region 
of our Galaxy. We study metallicities of 33 K/M giants situated in the Galactic Center region from observations obtained with the 
APOGEE survey. We selected K/M giants with reliable stellar parameters from the APOGEE/ASPCAP pipeline. Distances, interstellar 
extinction values, and radial velocities were checked to confirm that these stars are indeed situated in the inner Galactic Bulge. We 
find a metal-rich population centered at [M/H] = -1-0.4 dex, in agreement with earlier studies of other bulge regions, but also a peak at 
low metallicity around [M/H] = -1.0 dex, suggesting the presence of a metal-poor population which has not previously been detected 
in the central region. Our results indicate a dominant metal-rich population with a metal-poor component that is enhanced in the 
a-elements. This metal-poor population may be associated with the classical bulge and a fast formation scenario. 
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1. Introduction 

The Milky Way bulge is such a complex system that its for¬ 
mation and evolution are still poorly understood. Due to the 
high extinction and crowding, the study of the Galactic bulge 
remains challenging. Extinction of more than 30 mag in Ay 
in the Galactic Center (GC) regions requires IR spectroscopy. 
While more and more detailed chemical abundances in the in¬ 
termediate and outer Bulge (such as Baade’s Window) are now 
available thanks to the large spectroscopic surveys ARGOS 
(Freeman et al. 2013|), Gaia-ESO (Rojas-Arriagada et al.|2014| l. 


and APOGEE ( Garcia Perez et al.|2013| l, chemical abundances 

of stars in the Inner Galactic Bulge (IGB) with projected dis¬ 
tances of Rg < 200 pc from the GC remain poorly studied. 

Most of the chemical abundance studies in the IGB have 
been limited to luminous supergiants whose complex and ex¬ 
tended stellar atmospheres make abundance analysis difficult. 
|Carr et aTT ( 2000|l, [Ramirez et al.| ( |20(j0l l, |Cunha et al.| ( |2007| l, 
and Davies et al.j ( 2009| l analysed high-resolution spectra of 
supergiant stars in the GC region and found metallicities of 
near solar metallicity. A similar conclusion was obtained by 
[Najarro et al.j (|2009|l after analysing two luminous blue vari¬ 


ables. Cunha et al. ( 2007| l derived abundances nine giant and 
supergiant stars in the Central Cluster and found a metallicity 
of [Fe/H] = -I-0.14 + 0.06, together with enhanced [O/Fe] and 
[Ca/Fe] abundances. The total range in Fe abundance among GC 
stars, 0.16 dex, is significantly narrower than the iron abundance 
distributions found in the literature for the older bulge popula¬ 
tion. Ryde & Schultheis (2015, RS2015) studied nine field stars 
with a projected distance R < 50pc. They found a metal-rich 
population with [Fe/H] = -i-O.ll + 0.15 and low a-value, and a 
lack of a metal-poor population, similar to Cunha et al.j (2007 1 . 
Their mean metallicities are ~ 0.3 dex higher than fields in the 
inner bulge ( Rich et ar]|2007| [Rich et aL 2012|l, indicatin g that 
the GC region contains a distinct population. Grieco et al.| ( 2015] l 
compared these data with a chemical evolution model and con¬ 
cluded that in order to reproduce the observed [a/Fe] ratios, the 
GC region should have experienced a main strong burst of star 
formation and evolved very quickly with an IMF which contains 
more massive stars. All these studies comprise a limited num¬ 
ber of targets with precise chemical abundances. In this letter, 
we will discuss the metallicity distribution of 33 APOGEE K/M 
giants in the so-called “GALCEN field”, increasing significantly 
the statistics of K/M giant stars with reliable metallicities. 
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2. The Sample 

2.1. APOGEE 

The Apache Point Observatory Galactic Evolution Experiment 
(APOGEE; Majewski] 2012[ 2015 in prep.), one of four Sloan 
Digital Sky Survey-Ill (SDSS-III, Eisenstein et al.|2011| l exper¬ 
iments, is a large scale, near-IR, high-resolution (R ~ 22,500) 
spectroscopic survey of Milky Way stellar populations. The sur¬ 
vey uses a dedicated, 300-fiber, cryogenic spectrograph that is 
coupled to the wide-held, Sloan 2.5 m telescope ( |Gunn et al.| 
2006| l at the Apache Point Observatory (APO). APOGEE ob¬ 
serves in the //-band, where extinction by dust is signihcantly 
lower than at optical wavelengths (e.g., A(H)IA{V) ~ 0.16). 

With its high resolution and high S/N (~ 100 per Nyquist- 
sampled pixel), APOGEE determined both accurate radial ve¬ 
locities (better than 0.5 km s ') and reliable abundance measure¬ 
ments, including the most abundant metals in the universe (C, 
N, O), along with other a, odd-Z, and iron-peak elements. The 
latest SDSS-III Data Release (DR12; [Ahn et al.|2014[ ) provides 
spectra of more than 140,000 stars to the scientihc community, 
as well as the derived stellar properties, including radial veloc¬ 
ities, effective temperatures, surface gravities, metallicities, and 
individual abundances. Additional information, such as photom¬ 
etry and target selection criteria, is also provided and described 
in |Zasowski et al.] ( |2013| l. 


2.2. The GALCEN field 

The GALCEN APOGEE held is centered at (1, b)= (0.173°, 
-0.07°) with a held of view of ~ 2.9 sq. degree (see Eig. [^1, and 
it consists of 339 stars. Due to the high interstellar extinction 
(see, e.g., Schultheis et al.|2014 l, a special target selection pro¬ 
cedure has been applied for the GALCEN held (Zasowski et al. 
2013| l. Eorty eight stars are known to be Long-Period Variables 
based on K-band lightcurves (Matsunaga et al. 2009 1 . Thirty 
four stars are Asymptotic Giant Branch (AGB) stars based on 
low-resolution HK-band spectra (Schultheis et al. 2003| l. In addi¬ 
tion, spectroscopically known supergiants such as IRS7, IRS 19, 
and IRS22 were observed as well as supergiant candidates based 
on photometric criteria. These AGB and supergiant targets have 
among the coolest effective temperatures (Tgff < 3500K) in 
the APOGEE stellar sample. Eigure shows the distribution of 
the APOGEE GALCEN sample in Galactic coordinates. Only 
~ 10% of the 339 GALCEN stars (hlled circles) were used for 
our analysis (Sect. 3). 

Eigure shows the dereddened colour-magnitude diagram 
((H - Ko) vs. Kq) of the selected GALCEN sources (black 
points), together with the 2MASS data in this held that were 
dereddened using the high-resolution VVV extinction map of 
Schultheis et al. ( 2014| l. The GALCEN sample sources on the 
other hand, were dereddened using the RJCE method on a 
star-by-star basis (see [Zasowski et 2013 for more details). 
Superimposed are the Padova isochrones with an age of lOGyr 
and [Ee/H] = -0.4 and [Ee/H] = h-0.3. 


2.3. Stellar parameters 

Stellar parameters and chemical abundances of up to 15 ele¬ 
ments are determined by the ASPCAP pipeline. These values 
are based on a ;^^^-minimization between observed and synthetic 
model spectra (see Zamora et al.|2015 Holtzman et al.|20T5] for 
more details on the DR 12 spectral libraries) performed with the 
EERRE code (|Allende Prieto et al.||2006] and subsequent up- 
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Fig. 1. Galactic longitude vs. Galactic latitude distribution of 
the GALCEN sample superimposed on the 8 pm IRAC image 
of GLIMPSE-II ( jChurchwell eraLl|2009] l. The red hlled sym¬ 
bols indicate our selected sample with [M/H] < -0.7, the blue 
ones with [M/H] < -0.7 and the open black symbols the full 
GALCEN sample. Circles denote galactocentric radii of 50, 100, 
and 150pc, assuming Rq = 8kpc. The black asterisk show 
the GC held of Cunha et al. (2007), the black triangles those 
of Ryde&Schultheis (2015), black squares those of Rich et al. 
(2007). 


dates). Model spectra are interpolated on a regular grid com¬ 
puted with the ASSeT code ( [Koesterl^ 2009|l for a custom line 
list specially compiled for the surv ey ([Shetrone et al.|2015| l, and 
a large grid of model atmospheres ( M6szaros et al.|2()T2 l. 

The accuracies of the ASPCAP logg, and [M/H] val¬ 
ues were evaluated by Holtzman et al. ( 2015| l using a sample of 
well-studied held and cluster stars, including a large number of 
stars with asteroseismic stellar parameters from NASA’s Kepler 
mission ( jBorucki et al.|20f0| l. In this work, we use DR12 cali¬ 
brated stellar parameters (we refer to Holtzman et al.||20T5 for 
details). 


3. Sample selection 

In order to ensure the most reliable stellar parameters and abun¬ 
dances, we selected stars in the GALCEN held from the DR12 
dataset with following criteria: (i) ASPCAP < 20; (ii) 
PARAMPLAG=0 for Teff, logg, and [M/H], meaning that we 
avoid stars which are too close to the edge of the model grid; and 
(iii) log g < 3 to avoid foreground dwarfs. We exclude known su¬ 
pergiants, AGBs or possible AGB/supergiant candidates (based 
on near-and mid-IR criteria) for our analysis but restrict our¬ 
selves to K/M giants. In total we are left with 33 K/M giants. 
We checked by visual inspection the reliabilty of the ASPCAP 
hts to the observed spectra for each star. Due to the cool ef¬ 
fective temperatures of our stars (see Pig. [^, we restricted our¬ 
selves to using temperatures, surface gravities, metallicities, and 
global cr/Pe abundances, rather than the individual elemental 
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Fig. 2. Dereddened (H - K)o vs. Kq 2MASS diagram using the 
extinction map of Schultheis et al. (2014). The blue line shows 
the Padova isochrones of lOGyr with [Fe/H] = -0.4 while 
the green line [Fe/H] = -1-0.3. Black circles show our sample 
of stars dereddened with the RJCE method. Open circles show 
[Fe/H] > —0.7 while the filled ones are the low metallicity stars 
with [Fe/H] < -0.7. 


abundances. We added an additional cut in effective tempera¬ 
ture (Teff > 3700 K) for the analysis of the cr-elements to avoid 
possible systematic abundance effects in cool stars. 

Figure shows the distribution of the sample stars in Teff 
vs. logg space, colored by metallicity. One sees clearly that our 
sample covers a wide metallicity range and that the majority 
of our targets are M giants. As we exclude cooler stars with 
Teff < 3600 K, we could be missing more metal-rich stars. In 
order to be sure that our stars are really located in the 1GB, spec- 
trophotometric distances were calculated using the Bayesian ap¬ 
proach of Anders et al.| ( |2014) l. Our stars are indeed located at a 
heliocentric distance of 8 kpc with a typical uncertainty in the or¬ 
der of 20%. No difference in the distance distribution was found 
between metal-rich and metal-poor stars. In addition, our derived 
extinction values are compatible within 20% with those of the 
IGB. We also studied the radial velocity distribution of our sam¬ 
ple, which has a mean value of -32kms * with a dispersion 
of 122kms '. These values ar e similar to the v elocity disper¬ 
sion of 114kms“* obtained by 


Zhu et al. 


( 2008|l for the central 
parsec of the Gal axy and to the results obtained in the study of 


Rich et al.| ( |2007[ ), with a velocity dispersion of 115 + 15 km s“*. 
This is in perfect agreement with the predicted dispersion by the 
N-body bulge model of Shen et al. ( |2010j l. Interestingly, there 
is one metal-poor [Fe/H] = -1.0 high radial velocity star with 
Vrad = 356 km s '. Proper motions would be necessary to derive 
the orbital period of this star. 


4. Results and Discussion 

Figure]^ shows the histogram of the GALCEN sample metallic- 
ities. Eor comparison we show in green the sample of [Rich et al.| 


Fig. 3. Teff vs. logg as function of [M/H] for our selected sample 


( 2007| l and in blue, that of |Ryde & Schultheis ( 2015| l. We also 
show in red the metallicity distribution from the central cluster 
sample of Cunha et al. ( |2007| l, although this should be consid¬ 
ered a different stellar population. We clearly see a prominent 
peak in the metallicity distribution at ~ 0.4 dex, which is about 
0.2-0.3 dex more metal-rich than found in [Ryde & Schultheis^ 


( 2015| l and Rich et al. ( 2007| l. The most striking feature in the 
metallicity distribution is the large dispersion, with a r.m.s. scat¬ 
ter of 0.55 dex, which is much larger than seen in previous 
studies (~ 0.15 dex), and the presence of a metal-poor peak 
at [Fe/H] ~ -1.0, which has been not observed until so far 
in the IGB. This metal-poor population could be indicative of 
a classical bulge and a ssociated with a fast formation scenario 
( Babusiaux et al.|2010 i. A manual analysis of the '^C/'^C ratio 
in one of the metal-poor stars (2M17480557-2918376) reveals a 
value of '^C/'^C =5+2. Such a low carbon isotope ratio for a red 
giant star is indicative of a low-mass first ascent giant (M~0.8 


1.0 Mq) as first quantified by Gilroy 


1989 1 for open cluster 


giants and confirmed by Mikolaitis et al. (2012 1 . This value is 
close to that found in typical globular cluster red giants (e.g.. 


Briley et al. 

11994)1 which typify an old, low-mass population. 

Contrary to 

Ryde & Schultheis P015i, |Rich et al. P0071l, and 

Cunha et al. 

( 2007|), the metallicity distribution resembles a typ- 


etar] ( |20TT] i, etc. for bulge fields at higher latitudes (b > 4°). It 


should be noted that the metal-rich peak of our studied sample 
agrees with the metal-rich peak of the galactic bulge sample of 


Bensby et al. (20131 and the red clump bulge sample of Hill et al. 
( 201 Ij l. However, the second peak in their metallicity distribu¬ 
tion is at-0.5 to -0.6 dex (compared to [M/H]-1.0 in our 

sample) with a lack of metal poor stars with [M/H] < -1. Our 


derived metallicity dispersion is comparable to that of Bensby 
et al. (2013|l (0.53 dex) but larger than the red clump distribution 


of |Hili et al.| ( |201 1) 1 in Baade’s Window. 

What is the reason that this metal-poor population has not 
been revealed by other studies? The Ryde & Schultheis ( |2015| l 
and [Rich et aT] (|2007|l work contain a much smaller sample. 
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[M/H] 



Rich et al. 


Fig. 4. Histogram of the metallicities comp ared to _ 

(2007 I (green), Cunha et al. ( 2007| l (red) and jRyde & Schultheis 
(2015 I (blue). The bin size is 0.1 dex. 


Fig. 5. a vs [M/H]. The blue filled circles indicate the galactic 
bulge data from Gonzalez et al. ( 201 l| l, black from Bensby et al. 
(2013), green from Ryde&Schultheis (2015), and in red those of 
our GALCEN sources. We used only stars with Tgff > 3700 K. 


which could be one reason that they have just missed these 
metal-poor stars. They deredden their targets using global ex¬ 
tinction maps such as Schultheis et al.|('l999[), Schultheis et al. 


( |2009| l, or Gonzalez et al.| ( 2012| l. On the contrary, APOGEE 


calculates the reddening values on a star-by-star basis using a 
combination of near and mid-lR photometry and applies a ho¬ 
mogenous (J - K)o > 0.5 colour cut (see Zasowski et al.|[20T3 
for further details). 

Eigure [^displays the a-abundances of our sample stars as a 
function of [M/H]. As a comparison we also show the Galactic 
Bulge sample of the dwarfs and subgiants stars of|Bensby et al. 


(|2013|l, the M giants sample in the Galactic Center ofjRyde & 


Schultheis (|2015 1 , and the red clump star sample in Baade’ 


Window of|Gonzalez et al. (|201 Ij). O ur a-abundances agree with 
those of Ryde & Schultheis ( 2015| l, Gonzalez et al. pOll i and 
Bensby et al. j2013|l on the metal-rich side. However, we notice 


that the a-elements of |Bensby et al.| ( [2013| ) and [Gonzalez et aL| 
(|201 If are systematically higher than ours for low-metallicity 
stars ([Ee/H] < -0.5). 

None of the current models is able to reproduce the observed 
metal-poor peak seen in our metallicity distribution. [Grieco et al. 


(|2015|l calculated for the first time a chemical evolution model 


for the inner 200 pc and compared it with the dataset of Ryde 
& Schultheis (2015|l. In order to produce their observed [a/Fe] 


ratios and their metallicity distribution, they required a model as¬ 
suming a main early strong burst of star formation with a star for¬ 
mation efficiency of 25 Gyr“^ and a heavy IMF. However, their 
MDF consists of a single metal-rich population. 

Regarding the explanation for other metal-poor bulge 
peaks/populations, [Grieco et al. ( |2012| l, for example, used a 
chemical evolution model with two distinct populations that re¬ 
produces the MDF of [Hill et ah] ( [2011[ ) in Baade’s window. 
They concluded that the metal-poor population (which is dif¬ 
ferent than the one discovered here) was formed on a very short 
timescale associated with an intense burst of star formation and 


a high star formation efficiency, while the metal-rich population 
was formed on a longer time-scale and is related to the Galactic 
bar. |Tsujimoto & Bekki ( 2012| l modeled the MDF with a two- 
component disk model; a metal-poor component with a very 
short timescale of 1 Gyr and identical to the solar-neighbourhood 
thick disc, and a metal-rich component with a longer timescale 
and a top-heavy IMF. 

This shows that chemical evolution models can reproduce 
observed MDFs but fail to predict MDFs in different parts of 
our Galaxy. Thus, detailed confrontation with chemical evolu¬ 
tion models are necessary to understand the origin of the metal- 
poor peak. 
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